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Abstract 


The galactic heavy ion transport code (GCRTRN) and the nu- 
cleon transport code (BRYNTRN) are integrated into a code package 
(HZETRN). The code package is computer efficient and capable of op- 
erating in an engineering design environment for manned deep space 
mission studies. The nuclear data set used by the code is discussed 
including current limitations . Although the heavy ion nuclear cross 
sections are assumed constant , the nucleon-nuclear cross sections of 
BRYNTRN with full energy dependence are used. The relation of the 
final code to the Boltzmann equation is discussed in the context of sim- 
plifying assumptions. Error generation and propagation is discussed, 
and comparison is made with simplified analytic solutions to test nu- 
merical accuracy of the final result. A brief discussion of biological 
issues and their impact on fundamental developments in shielding tech- 
nology is given. 

1. Introduction 

Propagation of galactic ions through matter has been studied for the past 40 years as a means 
of determining the origin of these ions. Peters (1958) used the one-dimensional equilibrium 
solution and ignored ionization energy loss and radioactive decay to show that the light ions 
have their origin in the breakup of heavy particles. Davis (1960) showed that one-dimensional 
propagation is simplistic and that leakage at the galactic boundary must be taken into account. 
Ginzburg and Syrovatskii (1964) argued that the leakage can be approximated as a superposition 
of nonequilibrium one-dimensional solutions. The “solution” to the steady-state equations is 
given as a Volterra equation by Gloeckler and Jokipii (1969), which is solved to first order in 
the fragmentation cross sections by ignoring energy loss. They provide an approximation to the 
first-order solution with ionization energy loss included that is only valid at relativistic energies. 
Lezniak (1979) gives an overview of cosmic ray propagation and derives a Volterra equation 
including the ionization energy loss which he refers to as a solution “only in the iterative sense” 
and evaluates only the unperturbed term. No attempt is made to evaluate either the first-order 
or higher order perturbation terms. The main interest among cosmic ray physicists has been in 
first-order solutions in the fragmentation cross sections, since path lengths in interstellar space 
are on the order of 3 4 g/cm 2 . Clearly, higher order terms cannot be ignored in accelerator 
or space shielding transport problems (Wilson 1977a, 1977b, and 1983; Wilson et al. 1984). 
Aside from this simplification, the cosmic ray studies that were discussed have neglected the 
complicated three-dimensional nature of the fragmentation process. 

Several approaches to the solution of high-energy heavy ion propagation including the 
ionization energy loss have been developed (Wilson 1977a, 1977b, and 1983; Wilson et al. 
1984, 1989a, and 1987b; Wilson and Badavi 1986; Wilson and Townsend 1988; Curtis, Doherty, 
and Wilkinson 1969; Allkofer and Heinrich 1974; Chatterjee, Tobias, and Lyman 1976; Letaw, 
Tsao, and Silberberg 1983; Ganapol, Townsend, and Wilson 1989; Townsend, Ganapol, and 
Wilson 1989) over the last 20 years. All but one (Wilson 1977a) have assumed the straight 
ahead approximation and velocity conserving fragmentation interactions. Only two (Wilson 
1977a; Wilson et al. 1984) have incorporated energy- dependent nuclear cross sections. The 
approach by Curtis, Doherty, and Wilkinson (1969) for a primary ion beam represented the 
first generation secondary fragments as a quadrature over the collision density of the primary 
beam. Allkofer and Heinrich (1974) used an energy multigroup method in which an energy- 
independent fragmentation transport approximation was applied within each energy group after 
which the energy group boundaries were moved according to continuous slowing down theory 
(-dE/dx). Chatterjee, Tobias, and Lyman (1976) solved the energy- independent fragment 



transport equation with primary collision density as a source and neglected higher order 
fragmentation. The primary source term extended only to the primary ion range from the 
boundary. The energy-independent transport solution was modified to account for the finite 
range of the secondary fragment ions. 

Wilson (1977b) derived an expression for the ion transport problem to first order (first 
collision term) and gave an analytic solution for the depth-dose relation. The more common 
approximations used in solving the heavy ion transport problem were further examined by 
Wilson (1977a). The effects of conservation of velocity on fragmentation and the straight ahead 
approximation are found to be negligible for cosmic ray applications. Solution methods for 
representing the energy-dependent nuclear cross sections are developed (Wilson 1977a). Letaw, 
Tsao, and Silberberg (1983) approximate the energy loss term and ion spectra by simple forms 
for which energy- derivatives are evaluated explicitly (even if approximately). The resulting 
ordinary differential equations in position are solved analytically similar to the method of 
Allkofcr and Heinrich (1974). This approximation results in a decoupling of motion in space 
and a change in energy. In Letaw’s formalism, the energy shift is replaced by an effective 
attenuation factor. Wilson (1983) adds the next higher order (second collision) term. This 
term was found to be very important in describing 20 Ne beams at 670 MeV/amu. The three- 
term expansion was modified to include the effects of energy variation of the nuclear cross 
sections (Wilson et al. 1984). The integral form of the transport equation was also used to 
derive a numerical marching procedure to solve the cosmic ray transport problem (Wilson 
1977a; Wilson and Badavi 1986). This method can easily include the energy-dependent nuclear 
cross sections w-ithin the numerical procedure. Comparison of the numerical procedure with 
an analytical solution to a simplified problem (W ilson and Badavi 1986; W ilson and Townsend 
1988) validates the solution technique to about 1 percent accuracy. Several solution techniques 
and analytic methods have been developed for testing future numerical solutions to the transport 
equation (Ganapol, Townsend, and Wilson 1989; Townsend, Ganapol, and W ilson 1989). More 
recently, an analytic solution for the laboratory ion beam transport problem has been derived 
with a straight ahead approximation, velocity- conservation at the interaction site, and energy-- 
independent nuclear cross sections (W r ilson et al. 1989a). 

In this overview of past developments, the applications split into two separate categories 
according to a single ion species with a single energy at the boundary- versus a broad 
host of elemental types with a broad, continuous energy spectrum. Techniques requiring a 
representation of the spectrum over an array of energy- values require vast computer storage and 
computation speed for the laboratory beam problem to maintain sufficient energy- resolution. 
On the other hand, analytic methods (Wilson 1977a) are probably best applied in a marching 
procedure (Wilson and Badavi 1986), which again has within it a similar energy- resolution 
problem. This is a serious limitation because we require a final High Charge and Energy (IIZE) 
Code for cosmic ray shielding that has been validated by laboratory experiments. There is 
some hope of having a single code which is able to be validated in the laboratory (W ilson et al. 
1984; Schimmerling et al. 1989; Schimmerling 1990). More recently a Green’s function has been 
derived which holds promise for a code which may be tested in the laboratory environment and 
applied to space radiation protection (Wilson et al. 1990). 

In the present paper, we start with the general Boltzmann equation and simplify by using 
the standard assumptions to derive the straight ahead equation in the continuous slowing 
dow-n approximation with the assumption that heavy projectile breakup conserves velocity. 
A numerical procedure is derived with the coupling of the heavy ions to the nucleon fields. 
Numerical stability and error propagation are discussed. Analytic benchmark solutions are 
used to test the numerical procedures. The nuclear data and nuclear models are presented. The 
general computational procedures and data base have been in use for GCR shielding studies 
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since about 1987. A discussion of related biological issues and their impact on the need for 
future code development is also given. 

2. Boltzmann Equation 

Since the volume of any material is mostly electrons, it is natural that most of the interactions 
of energetic ions passing through any material are with electrons. The cross section for the 
interactions of atomic electron cr at is 


rr at - 1(T 16 cm 2 

The long range of the nuclear coulomb field also presents a sizable cross section a r to the passing 
ion 

<t c ~ 1(T 19 cm 2 

and ion collisions are dominated by these two processes but individual collisions have little effect 
on the passing ion. 

Although most collisions in the material are coulomb collisions with orbital electrons and 
nuclei, the rare nuclear reactions are of importance because of the large energy transferred in the 
reaction and the generation of new energetic particles. This process of transferring kinetic energy 
into new r secondary radiations occurs through several different processes, such as direct knockout 
of nuclear constituents, resonant excitation followed by particle emission, pair production, and 
possible coherent effects within the nucleus. Through these processes, a single particle incident 
on the shield may attenuate through energy transfer to electrons of the media or generate a 
multitude of secondaries causing an increase in exposure (transition effect). Which process 
dominates depends on energy, particle type, and material composition. This development 
of cascading particles is depicted in figure 2.1 as a relative comparison between high-energy 
proton and alpha particle cascades in the Earth’s atmosphere. Note the similarities displayed 
in figure 2.1 for individual reaction events and the nuclear star events shown in figure 2.2 for 
nuclear emulsion (Krebs 1950). 


Protons 




Figure 2.1. Cascade development in low Z materials. 
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(a) Z « 20 projectile, low-energy transfer event. 



• v ',y js 

(b) Neutral primary, 27 relativistic charged 
prongs. 



(c) Z ^ 17 projectile, high-energy transfer 
event. 


Figure 2.2. Nuclear star events observed in nuclear emulsion (from Krebs 1950). 



Figure 2.3. Transport of particles through spherical region. 




The relevant transport equations are derived on the basis of conservation principles. Consider 
a region of space filled by matter described by appropriate atomic and nuclear cross sections. In 
figure 2.3, we show a small portion of the region enclosed by a sphere of radius 6. The number 
of particles of type j leaving a surface element 8 2 dfl is given as (f>j(x + 8fl,fl ) E)8 2 dfl where 
<t>j(x } ft, E ) is the particle flux density, x is a vector to the center of the sphere, ft is normal to 
the surface element, and E is the particle energy. The projection of the surface element through 
the sphere center to the opposite side of the sphere defines a flux tube through which pass a 
number of particles of type j given as — 8fl,ft, E)6 2 dft , which would equal the number 
leaving the opposite face if the tube defined by the projection were a vacuum. The two numbers 
of particles in fact differ by the gains and the losses created by atomic and nuclear collisions as 
follows: 


4>j(x -f <5f2, ft , E)6 2 dft = (f>j(x — 6ft , ft , E)6 2 dfl 


+ <5 2 dft J S diY, j n ; , E, E') ^(x + m.O',#') dfl' dE' 

r6 

— 6 2 dft / diffj(E) <f>j(x + eft,ft,E) (2.1) 

J —6 


where crj(E) and aji i (ft,ft' , E, E') are the media macroscopic cross sections. The cross section 
E,E ! ) represents all those processes by which type k particles moving in direction 
ft' with energy E' produce a type j particle in direction SI with energy E . Note, there may 
be several reactions which may accomplish this result and the appropriate cross sections of 
equation (2.1) are the inclusive ones. Note that the second term on the right-hand side of 
equation (2.1) is the source of secondary particles integrated over the total volume 28{8 2 dft ) 
and the third term is the loss through nuclear reaction integrated over the same volume. We 
expand the terms of each side and retain terms to order <5 3 explicitly as 

6 2 dft [<t>j{x,fl,E) + 5n-V^j(x,n f JB)] = 6 2 dft [^(x, E) - 6fl ■ V^(x, ft, E) 

+ 2 6 W a jk {fl,ft f ,E,E f ) 4> k {x,fl\Ef) dfl f dE f 
k J 

-2 6<Tj(E) ^(x,n,E)] +0(<5 4 ) (2.2) 

which may be divided by the cylindrical volume 26(<5 2 dfl) and written as 

n • V^(x, fl,E) = Y, [ o', E, E')<j> k (x , ft E') dfl' dE' 

k J 

-aj(E) <j> j (x,n i E) + 0(8) (2.3) 

for which the last term 0(<5) approaches zero in the limit as 8 — ► 0. Equation (2.3) is recognized 
as a time-independent form of the Boltzmann equation for a tenuous gas. Atomic collisions 
(i.e., collisions with atomic electrons) preserve the identity of the particle, and two terms of the 
right-hand side of equation (2.3) contribute. The differential cross sections have the approximate 
form for atomic processes 

<r%fl, ft', E, E') = Y, - !) tjk 6{E + e n - E 1 ) (2.4) 

n 


where n labels the electronic excitation levels and e n represents the corresponding excitation 
energies which are small (1-100 eV in most cases) compared with the particle energy E. The 
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atomic terms may then be written as 

X / <rf k (n, ft', E, Ef) <t> k (x, ft', E') (in' dE' - crf(E) </>j(x, SI, E) 

k J 

= X a jri( E + f n) <^(x> Sl,E + f„) - (jf{E) SI, E) 
n 

« X <(£) ^( X - + X Hn( £ ) «. £0] 

n n 

- °f{E) 4>j{x,VL,E) 

= ^[Sj(E)<t>j(x,n,E)) 

since the stopping power is 

S j (E) = '£a jn (E)e n 

n 

and the total atomic cross section is 

af{E) - X<( E ) ( 2 ' ? ) 

71- 

Equations (2.5) to (2.7) allow us to rewrite equation (2.3) in the usual continuous slowing down 
approximation as 

n 

n • vd>j(x, n,E)~ — [Sj(e) 4>j(x, n, e)] + <tj(E) <^(x, si, e) 

= J X E > E ') M*, n', £') dn' dE' (2.8) 

where the cross sections of equation (2.8) now contain only the nuclear contributions. 

The purpose of the rest of this report concerns finding values for the atomic and nuclear cross 
sections, evaluating solutions to equation (2.8) for various boundary conditions, and making 
application to various radiation protection issues. 

The response of materials to ionizing radiation is related to the amount of local energy 
deposited and the manner in which that energy is deposited. The energy given up to nuclear 
emulsion is shown for several ions in figure 2.4. The figure registers developable crystals caused 
by the passage of the particle directly by ionization or indirectly by the ionization of secondary 
electrons (6-rays). These 6-rays appear as hairs emanating from the particle track. Note that 
the scale of the 6-ray track is on the order of biological cell dimensions (2 TO /mi). Many of the 
modern large integrated circuits are even of the 0.5 /xm scale. For this scale, track structure 
effects become important as interruptive events as a particle passes through active elements of 
such circuits. _ 

From the radiation protection perspective, the issues of shielding are somewhat clearly drawn. 
Given the complex external environment, the shield properties alter the internal environment 
within the spacecraft structure as shown in figure 2.5. The internal environment within the 
spacecraft interacts with onboard personnel or equipment. If sufficient knowledge is known 
about specific devices and biological response, then the shield properties can be altered so as 
to minimize adverse affects. Since the shield is intimately connected to the overall engineering 
systems and often impacts launch cost, the minimization of radiation risk is not independent 
of other risk factors and mission cost. Even mission objectives are at times impacted by radiation 


(2.5) 

( 2 . 6 ) 
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protection requirements (e.g., the Viking solar cell design to ensure sufficient solar power in the 
event of decreased performance caused by a large solar flare during the mission affected the 
weight allowed the experiments package). The uncertainty in shield specification is clearly an 
important factor when such critical issues are being addressed. There is uncertainty in subsystem 
response which can be easily (more or less) obtained for electronic or structural devices. The 
uncertainty in response of biological systems is complicated by the long delay times (up to 30 
years) before system response occurs and the unusually small signal-to-noise ratio in biological 
response. Clearly a difficult task remains before risk assignments can be made for long duration 
deep space missions. 

3. Transport Formalism 

The Boltzmann equation in the continuous slowing approximation is given by 




<t>j (x,n,E) 


]T j dE' dfl' a jk {E, E\ O, O') 4>k (*> &) 


(3.1) 


where <f>j (x, ft, E) is the flux of ions of type j with atomic mass Aj at x with motion along ft 
and energy E in units of MeV/amu (note that the energy units of eq. (2.8) are MeV), cr 3 {E) 
is the corresponding macroscopic cross section, Sj(E ) is the linear energy transfer (LET), and 
(F, F', ft, ft') is the production cross section for type j particles with energy E and direction 
ft by the collision of type k particle of energy E r and direction ft'. The term on the left-hand 
side of equation (3.1) containing Sj(E) is a result of the continuous slowing down approximation, 
whereas the remaining terms of equation (3.1) are the usual Boltzmann terms. The solutions 
to equation (3.1) are unique in any convex region for which the inbound flux of each particle 
type is specified everywhere on the bounding surface. If the boundary is given as the loci of the 
two-parameter vector function 7(5, t) for which a generic point on the boundary is given by T, 
then the boundary condition is specified by requiring the solution of equation (3.1) to meet 


<t> j (r,n,E) = ii> J (T,n,E) 


(3.2) 


for each value of ft such that 

ft-n(T) <0 (3.3) 

where n(T) is the outward-directed unit normal vector to the boundary surface at the point T 
and ipj is a specified boundary function. 

The fragmentation of the projectile and target nuclei is represented by the quantities 
&jk (F, F , ft, ft ) , which are composed of three functions: 

a jk (E, E\ O, O') = c k {E') v ]k {E') f jk (E, E' , O, O') (3.4) 

where ^(F') is the average number (which we loosely refer to as multiplicity) of type j particles 
being produced by a collision of type k of energy F' and f 3 ^ (F, F', ft,ft') is the probability 
density distribution for producing particles of type j of energy E into direction ft from the 
collision of a type k particle with energy E f moving in direction ft'. For an unpolarized source of 
projectiles and unpolarized targets, the energy angle distribution of reaction products may be 
a function of the energies and cosine of the production angle relative to the incident projectile 
direction. The secondary multiplicities i 'jk(E) and secondary energy angle distributions are the 
major unknowns in ion transport theory. 
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Information on the multiplicity v jk {E) was obtained in the past through experiments with 
galactic cosmic rays as an ion source, and the fragmentation of the ions on target nuclei was 
observed in nuclear emulsion (Cleghorn, Freier, and Waddington 1968). Such data are mainly 
limited by not knowing the identity of the initial or secondary ions precisely and by relatively 
low counting rates of each ion type. The heavy ion acceleration by machine makes a reduction 
in the uncertainty possible because large count rates can be obtained with known ion types. 
In addition, the accelerator experiments are providing information on the spectral distribution 
fjk {E, E ' , ft, f2 r ) which has not before been available (Heckman et al. 1972). 

The spectral distribution function is found to consist of two terms that describe the 
fragmentation of the projectile and the fragmentation of the struck nucleus as follows (Heckman 
1975; Raisbeck and Yiou 1975): 


a jk (E, E', ft, n') - a k (E') Vj k (E f ) fj k (E, E\ ft, ft') + vJ k (E') fj k [E, E f , ft, ft') 


(3.5) 


where v p k and f P k depend only weakly on the target and uj k and fj k depend only weakly on 

the projectile. Although the average secondary velocities associated with f p are nearly equal to 
the projectile velocity, the average velocities associated with f T are near zero. Experimentally, 
Heckman (1975) observed that 


fjk (e> e\ n, iV) 


- 

3/2 r 

m 

V2 E exp 

2n ( a fk) . 

- 

r -| 

3/2 r 

m 

/ \ 2 

V2E exp 

. 2 ^ i a ik) . 



(p-pO 2 

HO 2 


(nv'2m£;-nV2m£’ / ) 2 



where p and pi are the momenta per unit mass of j and k ions, respectively, and 


(3.6) 


fJ k {E,E\W) 


- 

3/2 

- 

- 

m 

E exp 


P 2 

— — , 

CM 

r 

CM 

J 

2 l 

B) 2 J 


(3.7) 


where a p k and aj k are related to the root mean square (rms) momentum spread of secondary 
products. These parameters depend only on the fragmenting nucleus. Feshbach and Huang 
(1973) suggested that the parameters a p k and aj k depend on the average square momentum of 
the nuclear fragments as allowed by Fermi motion. A precise formulation of these ideas in terms 
of a statistical model was obtained by Goldhaber (1974). 

4. Approximation Procedures 

4.1. Neglect of Target Fragmentation 


Using equations (3.5), (3.6), and (3.7) in the evaluation of the source term (j(x,fi,£) of 
equation (3.1) results in 


tj (x, n, E) = £ f dE' dW a k {E') <f) k (x, nf E') U(E') ff k (E, E\ Q, fT) 

k J 

+ vJ k fJ k (E,E\n,rt)} 

= c f (x,n,E) + Cj(x,n,E) (4.1) 
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where, as before, the superscripts P and T refer to fragmentation of the projectile and target, 
respectively. The target term is 


Cj (x,n,£)=X; 

k 


r 

3/2 

- 

m 

V 2 E exp 

- m.E 

CM 

r 

CM 

1 

tan 



dE' vJ k (E r ) a k (E') <f> k (x,ft',£') 


which is negligibly small for 



m 


Thus, for calculating the flux at high energy, 


(4.2) 


(4.3) 


^ (x,n ,e) « cf (x,n,^) 


(4.4) 


4.2. Space Radiations 

A convenient property of space radiations is that they are nearly isotropic. This fact, coupled 
with the forward peaked spectral distribution, leads to substantial reductions in the source term 
as follows: 


C f (x, ft,!?) f dE' dft' <J k (E') Vj k (E') 

k J 


-,3/2 



V2E' 


x exp 


(n\/2 m E - Els/lmE^ 


21 



d> k (x, ft', e') 


(4.5) 


If 4> k (x, ft', E') is assumed to be a slowly varying function of ft', one may seek an expansion 
about the sharply peaked maximum of the exponential function. Such an expansion is made by 
letting 

ft' = ft + (cos 9 — l)f2 -}- e^sin# (4.6) 


where 

cos 6 = ft * fl f (4-7) 


and 


q x n f 

e<t> = |ft X ft' 


(4.8) 


The flux may be expanded as 


d> k {x,n',E') =4> k (x,ft,E') + 




[(cos 6 - l)ft + sin 0] + 


(4.9) 
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Substituting equation (4.9) into equation (4.5) and simplifying result in 


Cf(x,n,E) ^JdE’ a k (E') u^(E') 


rn 


2n 




3/2 


j/2_ 

\[E' 


exp 


[\/2mE — \J 2rnE / ^j 

>(-5)’ 


21 


X < 


<t>k (x, n, e') 


d 




WO 


2 1 


2m\pEE f 


+ ... 


(4.10) 


The leading term of equation (4.10) is clearly a good approximation to the source term whenever 

2mE fl- m <h (x,n,E') , . 

(^) 2 

The leading term is equivalent to assuming that secondary ions are produced only in the direction 
of motion of the primary ions. For space radiations which are nearly isotropic, relation (4.11) 
is easily met, and neglect of higher order terms in equation (4.10) results in the usual straight 
ahead approximation. If the radiation is highly anisotropic, then relation (4.11) is not likely to 
apply. Validity of the straight ahead approximation was discovered empirically by Alsmiller et al. 
(1965) and Alsmiller, Irving, and Moran (1968) for proton transport. 

4.3. Velocity Conserving Interaction 

Customarily, in cosmic ion transport studies (Curtis and Wilkinson 1972), the fragment 
velocities are assumed to be equal to the fragmenting ion velocity before collision. The order of 
approximation resulting from such an assumption is derived. Since the projectile energy E f is 
assumed to be equal to the secondary energy plus a positive quantity e, 


E f = E + e (4.12) 

and e is assumed to contribute to equation (4.10) only over a small range above zero energy, 
substituting equation (4.12) into equation (4.10) and expanding the integrand result in 


Cf (x,n,E) = yVfc(.E’) uf k {E) 

k 


(x,n,E) 


i 



nmE 


E— <p k (x, ft, E) 


\ 


(4Y 

nmE 


o 


2 mE 


+ 


(4.13) 


Because y \ a jk) i m ^ ^ 1 at those energies for which most nuclear reactions occur, the 

assumption of velocity conservation is clearly inferior to a straight ahead approximation but may 
be adequate for space radiations where the variation of ^(x,fl,F) with energy is sufficiently 
smooth. That is, 


E M 4>k(*,n,E) 


<Pk (x,ft,£) 
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4.4. Decoupling of Target and Projectile Flux 

Equation (3.1) with the use of equation (4.1) may be rewritten as 

(x, ft, E) = Fjk <f>k (X, ft, E) + F f $ k ( x > n > E ) 

k k 

where the differential operator is given by 


B J = 


« v - Tj ^ s AE) + °j(E) 


and the integral operator ^ Fj * = Fj k + Fj^j is given by 

F jk 4>k (x, n,E) = J dE' dfl' <r 3k (E, E\ ft, ft') <j> k (x, ft', Ef) 

Defining the flux as a sum of two terms 

4> 3 (x, ft, e) = 4>J (x, ft, E) + ( x , ft, e) 

allows the following separation: 

Bj 4>f (X, ft, E) = ]T Ff k 4$ (X, ft, E) + Y, Ffk $ (x, ft, E) 
k k 

Bj 4>J (x, ft, E) = ]T Ff k <$>1 (x, ft, E) + £ Ff 4% (x, ft, E) 

fc * 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.18) 

(4.19) 


As noted in connection with equations (4.1) through (4.4), the source term on the right-hand 
side of equation (4.19) is small at high energies and one may assume 

*J(x,ft,E)« 0 (4.20) 

for E (af) 2 /m. As a result of equation (4.20) and the fact that the ion range is small 
compared with its mean free path at low energy, one obtains 

Bj 4>f (x, ft, E) « Y, Ffk $ (x, ft, E) (4.21) 

k 

Bj 4>J (X, ft, E) « £ Ff (x, ft, E) (4.22) 

fc 

The advantage of this separation is that once equation (4.21) is solved by whatever means 
necessary, then equation (4.22) can be solved in closed form. The solution of equation (4.22) is 
accomplished by noting that the inwardly directed flux <pj must vanish on the boundary so that 


4>J (x, n,E)*£ dE' ^)sf(E) / dE " dQ (£,) £ " 
X 4>f (x + [Rj{E) - Rj(E')} ft, ft', E"} 
where E 7 = 1 [d + Rj(E)] with d the projected distance to the boundary. 


(4.23) 


12 


Using equations (3.5) and (3.7) in equation (4.23) yields 




(x,n,£) 


£ Pj(E) Sj(E) 


m 


2; r 


ten 


3/2 


n/ZE 7 


exp 


-m£' 


x<J'{T+[R j (E)-« J (£')]n} 


where 


Cj (*) = E / dE ' dn ' E ') t'l'k(E') 4>k (x,n',E') 


(4.24) 


(4.25) 


and crTJ. has been assumed to be a slowly varying function of projectile type k and projectile 
energy E. If the range of secondary type j ions is small compared with their mean free path 
lengths and the mean free paths of the fragmenting parent ions that is, 


Ri 


(-30 s 


77 } 


<4 


(4.26) 


then the integral of equation (4.24) may be simplified as 


4>j (x, fl, E) 


Cf(x) I*" 

Je 


Sj(E) 


777 


3/2 


(4) J 

which may be reduced into terms of known functions. Thus, 

3 m E 


VYi T' 


exp 


-m£' 

[Wf 


dE' (4.27) 


<t>i (x,n,E) 


^ #(*) 1 


Sj(E) 


2'K\pK 


2’ ter 


ten 


3 mE 


2 ’ ter 


te) J 


(4.28) 


in terms of the incomplete gamma function. One can show that equation (3.28) is equivalent to 




r 

— erfc 


mE 

erfc 


2 

A 

ten 

2 

A 


mE \ 


te) J 


T 


\ 



" - 



” 


mE 

—mE 


mE 1 

— mEcr 

\ 

( T \ 2 CXP 

( T \ 2 


( T \ 2 CXP 

( T \ 2 

? 

* w 

Lvw J 

M 

* bd 

Ivjk) \ 

J 


(4.29) 


At points sufficiently removed from the boundary such that 

,2 


J m 


(4.30) 
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equation (4.28) may be reduced to 


<t>j (x,ft 5 £) 


Sj(E) 


(*) 


2?r 


erfc 


mE 


+ 


mE 


N * Hi) 


exp 


—mE 


■tY 


K< 


(4.31) 


Jkj J 


The solution of equation (4.21) is now examined further. 


4.5. Back-Substitution and Perturbation Theory 

One approach to the solution of equation (4.21) results from the fact that the multiple- 
charged ions tend to be destroyed in nuclear reactions. Thus, 

F[ k = 0 (j > k) (4.32) 

This means that there is a maximum j such that 

Bj^(x,n,E) = 0 (4.33) 

where J is the largest j. Furthermore, 

Bj- 1 4> P J-X (x, n, E) = F?_ ltJ (X, n, E) (4.34) 

and, in general, 

N-\ 

Bj-n <t>j-N ( x > E) “ 5Z Fj-N,J-k fij-k ( x > B) (4.35) 

k~] 

for N < J - 1. Note that equations (4.34) and (4.35) constitute solvable problems. The singly 
charged ions satisfy 


B\ 4>f (x, a e) = Ff, 0 f (x, n, F) + F‘] k 4 (x, n, e) 


(4.36) 


k—2 


which, unlike equations (4.33) to (4.35), is an integral-differential equation that is difficult to 
solve directly. Equation (4.35) is solvable by perturbation theory, and the resultant series is 
knowm to converge rapidly for intermediate and low energies (Wilson and Lamkin 1975; Wilson 
et al. 1989a and 1989b; Wilson and Townsend 1988). Note that equations (4.33) and (4.35) are 
also obtained from perturbation theory as applied to equation (4.21) at the outset. Thus, the 
perturbation series is expected to converge after the first J plus a few terms. 


5. Galactic Ion Transport 

In the present section, wo expand on the methods developed earlier for nucleon transport 
(Wilson and Lamkin 1975) by combining analytic and numerical tools. The galactic cosmic 
ray ion transport problem is transformed to an integral along the characteristic curve of that 
particular ion. As a result of the conservation of velocity in fragmentation, the perturbation 
series (Wilson and Lamkin 1975) is replaced by a simple numerical procedure. The resulting 
method reduces the difficulty associated with the low-energy discretization and the restriction 
to a definite form for the stopping power. The resulting numerical computation is simple and 
nondemanding from computer requirements and yet gives superior results compared with other 
methods. 
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In the present work, we use the straight ahead approximation and neglect the target 
secondary fragments (Wilson 1977a and 1983). For multiply charged ions, the transport equation 
may be written as 


’ d_ 

dx 


dE 


Sj(E) + aj 


E) = J2 

k >j 


Tn jk a k 4>k{ X ,E) 


(5.1) 


where <pj(x ) E) is the flux of ions of type j with atomic mass Aj at x moving along the x axis at 
energy E in units of MeV/amu, a } is the corresponding macroscopic nuclear absorption cross 
section, Sj(E) is the change in E per unit distance, and m jk is the multiplicity of ion j produced 
in collision by ion k. The corresponding nucleon transport equation is 


0_ 

dx 


dE 


Si(E) + Oj (E) 


<t> 


<j(x, E) = Y, [ X a jk (E, E') 4> k (x, E') dE' 

l. JE 


(5.2) 


The quantities rrij ^ and crj are assumed energy independent in equation (5.1) but are fully energy 
dependent in equation (5.2). 

The range of the ion is given as 


f E 

Rj(E) = 

Jo 

The solution to equation (5.2) is found to be subject to boundary specification at x = 0 and 
arbitrary E as 

4> j {0 J E) = F j {E) (5.4) 

Usually, Ej(E) is called the incident beam spectrum. 

From Bethe’s theory, 

Sj(E) = S P (E) (5.5) 

and holds for all energies above 100 keV/amu provided the ions remain fully stripped. It follows 
that 

Z 1 Z v 

/ Rj(E) = ^Rp(E) (5.6) 

Aj A p 


dEJ_ 

Sj(E') 


(5.3) 


The subscript p refers to proton. Equation (5.6) is quite accurate at high energy and only 
approximately true at low energy because of electron capture by the ion which effectively reduces 
its charge, higher order Born corrections to Bethe’s theory, and nuclear stopping at the lowest 
energies. Herein, the parameter vj is defined as 

Vj Rj{E) = v k R k {E) (5.7) 


so that 



(5.8) 


Equations (5.6) to (5.8) are used in the subsequent development, and the energy variation in vj 
is neglected. The inverse function of Rj(E) is defined as 


E^RpiRjiE)} 


(5.9) 
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A method of solution is now discussed. For the purpose of solving equation (5.2), define the 
coordinates 

r]j = x — Rj(E) (5.10) 

tj = x + Rj{E) (5.11) 

where rfj varies along the particle path and is constant along the particle trajectory. The new 
fluence functions are taken as 


where 


Xji’Hj i £j) = Sj(E) 4>j{ x >E) — x Pj( x > r j) 

(5.12) 

Xk(Vj,£j) = Xk(Vk^k) 

(5.13) 

+ Vj — £>k + r lk 

(5.14) 

7j — £ j = ~ (Vk ~ £k) 

(5.15) 


and rj = Rj(E). Under this coordinate mapping, equation (5.1) becomes 


2 * b 

dm 


k k 


(5.16) 


where is assumed to be energy independent. There is a small variation in (»20 percent) 
which must eventually be taken into account. Solving equation (5.16) by using line integration 
with an integrating factor 

'1 


Hjivjitj) = ex P 


2°j«7 + Vj) 


(5.17) 


results in 


= exp 


~2 ff j($ + 7j) 


Xji £j’£j) 


i m 

+ - J exp 


2 J Hj 


2 Cr i( r ? Vj) 'y ] m jk a k XkiVk-^k^ ^T] 

Z J k k 


(5.18) 


where 


and 


j = V A + V U | V A - v ± c 

Vk 2u k n 2u k ^ 


r, V k ~ "j , ** + Vj e . 

2u k V + 2, k * 


ipj(x > r j ) = Xj(rij,Zj) 


Defining 
one may show 

i f j( x > r j) = vxp(-<rjx)i/>j(0,rj + x) + f dz exp(-ajz) J- ^ (x - z,r k + -j- 

J 0 h. k \ k 


(5.19) 


* (5.20) 
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Furthermore, it is easy to show that 


ipj(x + h, rj) — exp (-crjh) ipj(x,rj + h) 

rh . 

+ dz exp {-(JjZ) m jk a k ~ Tpk(x + h- z,r k + 

JO » 


u -i z ) 

V k ) 


(5.21) 


It is clear from equation (5.20) that 

ip k (x + h - z, r k ) = exp[-<7*.(/i - *)] ^ k (x, r k + h) + 0(h - z) 
which upon substitution into equation (5.21) yields 


(5.22) 


ipj{x + h, r j ) = exp(-ajh) ipj(x, rj + h) 

r h ^ . 

+ dz exp {-crjz)22 mjkak- 2 - exp [-a k (h - z 
JO k 


> 0 v k 

which is correct to order /i 2 . This expression may be further approximated by 
ipj{x + h, rj) = exp(-ajh) tpj(x, rj + h) 


)] 'Pk (z- T k + ~z + h- zj (5.23) 


+ E 


Vk 


exp (— (jj/i) — exp(— G^h) 
®k ~~ &j 


A (^, 


*■* + 


\s -I 

-*-h 

y k 


(5.24) 


which is accurate to O [(v k - Vj) h]. Equation (5.24) is the basis of the Galactic Cosmic Rays 
(GCR) Transport Code GCRTRN (Wilson and Badavi 1986; Wilson, Townsend, and Badavi 
1987a; Wilson and Townsend 1988). The nucleon transport equation (4.2) is solved by adding 
the heavy ion collision source of nucleons to the BRYNTRN code (Wilson et al. 1989). 

There are several quantities of interest that are now given. The integral fluence is given as 


The energy absorption per gram is 

Dj(x,>E) = f Ajipj[x, Rj(E)] dE (5.26) 

J E 

with the dose equivalent given as 

roc 

Hj{x,>E) = Aj / Q f iP 3 {x,R 3 {E)) dE (5.27) 

J E 

These quantities are used in shield design studies for protection against galactic cosmic rays. 

6. Analytic Benchmarks 

In this section, we address the question of GCR transport code validation. Ideally, validation 
should be accomplished with detailed transport data obtained from carefully planned and 
controlled experiments; unfortunately, there exists a paucity of such data. Although useful 
for comparison purposes, the atmospheric propagation measurements used previously (Wilson, 
Townsend, and Badavi 1987a) are clearly not definitive because they consist of integral fluences 
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of as many as 10 different nuclear species combined into a single datum. Although limited 
quantities of HZE dosimetry measurements from manned space missions (e.g., Skylab) are 
also available (Benton, Henke, and Peterson 1977), numerous assumptions concerning the 
relationships between dosimeter locations and spacecraft shield thicknesses and geometry must 
be made to estimate astronaut doses with GCR codes. Because many of these assumptions 
may involve inherently large uncertainties (a factor of 2 or greater), it becomes difficult to 
attribute any differences to particular assumptions or approximations that may have been used 
in the analyses. Without definitive GCR transport measurements with which to compare code 
predictions, other methods of validation must be considered. As noted by Wilson (1983) and 
Wilson et al. (1989a), there are several different versions of HZE transport codes available. 
When used with the same input spectra, interaction parameters, and boundary conditions, 
all should yield comparable results. The history of transport code development, however, 
suggests otherwise. For this reason, a realistic, nontrivial, exact, analytic solution to the 
simplified Boltzmann equation used to describe HZE transport has been formulated as an 
absolute standard for code comparison purposes. 

For the benchmark problem, the incident spectrum is limited to a single ion type ( j — </)• 
Because the GCR spectrum for a typical ion is of the form 

F(F)~£T 0 (6-1) 


where a « 2.5, we choose the energy spectrum to be of similar functional form as 


Fj(E) = 


OjJ 


[Rj(E)}' 


Sj{E) 


Defining the characteristic variables as 


and 


rij = x - Rj{E) 


£j = x + Rj{E ) 


( 6 . 2 ) 


(6.3) 


(6.4) 


equation (5.1) can be solved by the method of characteristics (Wilson 1977a; Wilson, et al. 
1989) to give 

*Ax,E)= . CXP( ~F.1... (6 5) 


where 

and 


[ujx + Fj(F)] 2 
E) = Sj(E) <Pj(x,E ) 


(6.6) 


aj = n/(l - mjj) ( 6 - 7 ) 

This is the trivial solution for the incident beam species. For j < J (secondary fragments), it 
can be shown that 


i’j(x,E) = ajrrijj ~ Ij(x,E) exp 


+ °J_ <;) 
2 


(6.8) 


where in terms of the exponential integral function E 2 (x) (see Abramowitz and Stegun 1964), 

exp \-b{uj + Vj)£j/2] 


Ij(x,E) 


VJ - Uj 


vfy vjtj 


(6.9) 
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for j = J — 1 and 


V J ~ v i 

Clearly, equations (6.9) and (6.10) are true for all j if rajy = 0 for all values of j < J (i.e., if the 
secondary fragments themselves do not fragment). 

The benchmark solution was calculated for an incident iron beam ( J = 26) in an aluminum 
target, for which the input parameters are <726 = 0.04568 cm 2 /g, <725 = 0.04260 cm 2 /g, and 
77725 26 026 = 0.00403 cm 2 /g. Results of the GCR transport code simulation of this benchmark 
for the propagating incident iron beam and secondary manganese (j = 25) ions and the exact 
analytic predictions obtained from equations (6.5) and (6.8) are given in tables 6.1 and 6.2. It is 
clear from these tabulated results that the numerical solution methods developed previously 
(Wilson and Badavi 1986; Wilson, Townsend, and Badavi 1987a) are accurate in solving 
equation (5.1) for GCR transport to within about 1 percent. This indicates that any limitations 
to accurately solving GCR transport problems must focus upon the simplifying approximations 
used to obtain equation (5.1) as well as upon unresolved issues concerning the need to include 
multiple coulomb scattering effects, fragment momentum dispersion effects, and perhaps most 
important, the nature and quality of the input cross-section data bases. To illustrate this point, 
we are aware of only one heavy ion transport code (Wilson et al. 1984), which uses energy- 
dependent cross sections. Recent studies, however, suggest that fully energy-dependent cross 
sections may be important for some transport code applications (Townsend and Wilson 1988). 


Table 6.1. Benchmark Numerical Simulation and Analytic Solution for Iron Ions 
as Function of Ion Depth" and Energy Into Aluminum Absorber 


E, 

MeV /amu 

t/>p e (0, E ) 

4'Fc( 10 > E ) 

^'FeC 20 ’ E ) 

Numerical 

Analytic 

Numerical 

Analytic 

Numerical 

Analytic 

0.0198 

0.1147 

1.090 

10.07 

100.1 

1059 

10 490 

1.394 X 10 s 
1.692 x 10 4 
9.217 x 10 2 
1.062 x 10 1 
9.310 x 10~ 3 
5.089 x 10“ 6 
2.970 x 10~ 8 

1.394 x 10 5 
1.692 x 10 4 
9.217 x 10 2 
1.062 x 10 1 
9.310 x 10 -3 
5.089 x 10 _G 
2.970 x 10~ 8 

4.334 X 10“ 5 
4.334 x 10“ s 
4.333 x 10” 5 
4.321 x 10“ 5 
3.699 x 10“ 5 
2.014 x 10 _G 
1.833 x 10“ 8 

4.382 x 10 -5 
4.381 x 10“ 5 
4.379 x 10 5 
4.360 x 10“ 5 
3.718 x 10“ 5 
2.019 x 10“ G 
1.833 x 10“ 8 

6.942 x 10“ G 
6.942 x 10“ G 
6.942 x 10 G 
6.932 x 10“ c 
6.400 x 10” 6 
8.741 x 10~ 7 
1.132 x 10“ 8 

7.044 x 10 _G 
7.044 x 10 _G 
7.043 x 10 _G 
7.027 x 10 _G 
6.478 x 10 _G 
8.799 x 10 -7 

1.132 x 10 -8 


"Depth is given in g/cm 2 . 


Table 6.2. Benchmark Numerical Simulation and Analytic Solution for Secondary Manganese 
Ions as Function of Ion Depth" and Energy Into Aluminum Absorber 


E, 

MeV /amu 

tW™. E) 

%n( 20 > 

Numerical 

Analytic 

Numerical 

Analytic 

0.0198 

0.1147 

1.090 

10.07 

100.1 

1059 

10490 

1.772 x 10~ G 
1.772 x 10 -6 
1.772 x 10“ 6 
1.767 x 10~ 6 
1.504 x 10 -6 
7.797 x 10 -8 
7.004 x 10- 10 

1.780 x 10~ G 
1.780 x 10 -6 
1.779 x 10“ 6 
1.771 x 10-® 
1.503 x 10~ 6 
7.806 x 10- 8 
7.004 x 10“ 10 

5.704 x lO -7 
5.704 x lO’ 7 
5.704 x 10- 7 
5.696 x 10~ 7 
5.242 x 10“ 7 
6.880 x lO' 8 
8.728 x 10“ 10 

5.768 x lO -7 
5.768 x lO -7 
5.767 x 10~ 7 
5.753 x 10- 7 
5.219 x 10~ 7 
6.918 x 10- 8 
8.728 x 10' 10 


"Depth is given in g/cm 2 . 
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7. High Charge and Energy (HZE) Nuclear Data Base 

The nuclear cross sections for neutron and proton interactions are described extensively by 
Wilson et al. (1989a and 1989b). The heavy ion absorption cross sections cr abs are currently 
derived from 

P , A t ) = nr 2 0 [4 /3 + A, - 0.4] 2 (7.1) 

which was fit to the asymptotic nuclear cross sections calculated by Wilson and Townsend 
(1981). 

In the abrasion-ablation fragmentation model, the projectile nuclei, moving at relativistic 
speeds, collide with stationary target nuclei. In the abrasion step, those portions of the nuclear 
volumes that overlap are sheared away by the collision. The remaining projectile piece, called a 
prefragment or primary residue, continues its trajectory with essentially its precollision velocity. 
As a result of the dynamics of the abrasion process, the prefragment is highly excited and 
subsequently decays by the emission of gamma radiation and/or nuclear particles. This step 
is the ablation stage. The resultant isotope, sometimes referred to as a secondary product, is 
the nuclear fragment whose cross section is measured. The abrasion process can be analyzed 
with classical geometric arguments (Bowman, Swiatecki, and Tsang 1973) or methods obtained 
from formal quantum scattering theory (Townsend et al. 1986a and 1986b). The ablation stage 
can be analyzed from geometric arguments (Bowman, Swiatecki, and Tsang 1973) or more 
sophisticated methods based upon Monte Carlo or intranuclear cascade techniques (Gosset 
et al. 1977; Hiifner, Schafer, and Schurmann 1975; Morrissey et al. 1978; Guthrie 1970). 
Predictions of fragmentation cross sections can also be made with the approximate semiempirical 
parameterization formulas of Silberberg, Tsao, and Shapiro (1976) and Silberberg, Tsao, and 
Letaw (1983). 


The amount of nuclear material stripped away in the collision of two nuclei of radius Rp and 
Rp is taken as the volume of overlap region times an average attenuation factor. The relevant 
formula for the constituents in the overlap volume in the projectile is given by the following 
formula; 


Aabr = FAp 




(7.2) 


where Cp and Cp are the maximum chord lengths of the intersecting surface in the projectile 
and the target, respectively, A is the nuclear mean free path, and the expressions for F differ 
depending on the nature of the collision (peripheral versus central) and the relative sizes of the 
colliding nuclei. 


For Rp > we have (Gosset et al. 1977) 


and 


P = 0.1250 w) 1 / 2 (i - 2 ) (^-y^) 2 


-0.125 


F = 0.75(1 - v) x ! 2 (t— -^) 2 

- 0.125 [3(1 - v) 1 ! 2 - l] (t^)‘ 


(7.3) 


(7.4) 
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with 


(7.5) 


V — 


P = 


Rp 

Rp + Rp 

b 

Rp + Rp 


and 


_ _ j _ R-t 

^ v Rp 


(7.6) 


(7.7) 


Equations (7.3) and (7.4) are valid when the collision is peripheral (i.e., the two nuclear volumes 
do not completely overlap). In this case, the impact parameter b is restricted such that 


Rp — Rp < b < Rp -K Rp 


(7.8) 


If the collision is central, then the projectile nucleus volume completely overlaps the target 
nucleus volume (b < Rp — Rp), and all the projectile nucleons are abraded. In this case, 
equations (7.3) and (7.4) are replaced by 


and 


P= -1 


F = 1 


(7.9) 

(7.10) 


and there is no ablation of the projectile because it was destroyed by the abrasion. 

For the case where Rp > Rp and the collision are peripheral, equations (7.3) and (7.4) 
become (Morrissey et al. 1978) 

-o. ™-*(H(W 

( V /2 (l _ 2 \ [( 1 / 0(1 - M 2 ) 1/2 - 1][(2 - mV ] 1/2 1 

kP J 


0.125< 0.5 ( - 




/ 


(- 


and 


F = 0.75(1 - i/) 1/2 


-& 


1-0 


— 0.125< 


3(1 - v) x R [1 - (1 - ^ 2 ) 3/2 ][l - (1 -mnfi-A 


(7.11) 


(7.12) 


where the impact parameter is restricted such that 

Rp — Rp < b < Rp + Rp (7.13) 

For a central collision (6 < Rp — Rp) with Rp > Rp , equations (7.11) and (7.12) become 




1/2 


and 


F = 1 - (1 


V> 3/2 ] [>-(!)' 


2i 1/2 


(7.14) 


(7.15) 


The charge ratio of removed nuclear matter is assumed to be that of the parent nucleus. 
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The surface distortion excitation energy of the projectile prefragment following abrasion of 
m nucleons is calculated from the clean-cut abrasion formalism of Bowman, Swiatecki, and 
Tsang (1973). For this model, the colliding nuclei are assumed to be uniform spheres of radii 
R t (i — P y T ). In collision, the overlapping volumes shear off so that the resultant projectile 
prefragment is a sphere with a cylindrical hole gouged out of it. The excitation energy is then 
determined by calculating the difference in surface area between the misshapen sphere and a 
perfect sphere of equal volume. This excess surface area AS is given by Gosset et al. (1977) as 

AS = 47r/?p[l + P — (1 — F) 2 / 3 ] (7.16) 

where the expressions for P and F differ depending upon the nature of the collision (peripheral 
versus central) and the relative sizes of the colliding nuclei which were given in previous 
equations. 

The excitation energy associated with surface energy is known to be 0.95 MeV/fm 2 for near 
equilibrium nuclei so that 

E f s = 0.95 AS (7.17) 

for small surface distortions. When large numbers of nucleons are removed in the abrasion 
process, equation (7.17) is expected to be an underestimate of the actual excitation. We therefore 
introduce an excess excitation factor in terms of the number of abraded nucleons A abr as 


IQ^abr , 2°^abr 
Ap A 2 p 


(7.18) 


which approaches 1 when the impact parameter is large but increases the excess excitation when 
many nucleons are removed in the collisions and when grossly misshapened nuclei are formed. 
The total excitation energy is then 

E s = Kf (7.19) 

which reduces to equation (7.17) for small A abr . We assume that all fragments with a mass 
of 5 are unbound, 90 percent of the fragments with a mass of 8 are unbound, and 50 percent of 
fragments with a mass of 9 ( 9 B) are unbound. 

A secondary contribution to the excitation energy is the transfer of kinetic energy of relative 
motion across the intersecting boundary of the two ions. The rate of energy loss of a nucleon 
when it passes through nuclear matter (Westfall et ah 1979) is taken at 13 MeV/fm, and the 
energy deposit is assumed to be symmetrically dispersed about the azimuth so that 6.5 MeV/fm 
per nucleon at the interface is the average rate of energy transfer into excitation energy. This 
energy is transferred in single particle collision processes, and the energy is transferred to 
excitation energy of the projectile for half of the events and leaves the projectile excitation 
energy unchanged for the remaining half of the events. The first estimate of this contribution is 
to use the length of the longest chord C\ in the projectile surface interface. This chord length 
is the maximum distance traveled by any target constituent through the projectile interior. 
The number of other target constituents in the interface region may be found by estimating 
the maximum chord C\ transverse to the projectile velocity which spans the projectile surface 
interface. The total excitation energy from excess surface and spectator interaction is then 


£i = 13Ci + ll3Ci(C t -1.5) 


(7.20) 


where the second term only contributes if Ct > 1.5 fm. We further assume that the effective 
longitudinal chord length for these remaining nucleons is one third the maximum chord length. 


i 

i 
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The decay of highly excited nuclear states is dominated by heavy particle emission. In the 
present model, we assume a nucleon is removed for every 10 MeV of excitation energy as 


A — ^ s 

abl “ 10 MeV 


(7.21) 


In accordance with the previously discussed directionality of the energy transfer, 
values as 


Ex — 




E x has two 


( 7 . 22 ) 


where Pj is the corresponding probability of occurrence of each value in collisions. 


The number of nucleons removed through the abrasion- ablation process is given as a function 
of impact parameter as 


A A = A abr (6) + A ab i(6) 


(7.23) 


The values of A A for carbon projectiles on a copper target and for copper projectiles on a carbon 
target are shown in figure 7.1. For each projectile, the dashed curve corresponds to E x = 0, 
whereas the solid curve corresponds to E x = E f x as given by equation (7.20). A real collision 
would be given by a statistical distribution between the limits shown by these two curves. The 
average event is calculated as if the two extremes occurred with equal probability, as noted in 
equations (7.22). 

The nuclear fragmentation parameters herein are approximated according to the abrasion- 
ablation model of Bowman, Swiatecki, and Tsang (1973). The cross section for removal of A A 
nucleons is estimated as 


cr(AA) = 7T&2 — 7T&i (7*24) 

where 62 is the impact parameter for which the volume of intersection of the projectile contains 
A abr nucleons and the resulting excitation energies release an additional A ab j nucleons at the 
rate of 1 nucleon for every 10 MeV of excitation such that 

Aabr(&2) + A abl (6 2 ) = A A- 1 - (7.25) 

and similarly for b\ 

A a br(&i) + A ab i(6i) = A A + - (7.26) 

The charge distributions of the final projectile fragments are strongly affected by nuclear 
stability. We expect that the Rudstam (1966) charge distribution for a given <r(AA) to be 
reasonably correct as 

a{A F , Z F ) = Fi exp (-R\Z F - SA F + TA 2 F | 3 / 2 ) <r(A A) (7.27) 

where R = 11.8 /A F , D = 0.45, S — 0.486, and T = 3.8 x 10 -4 according to Rudstam and Fj is 
a normalizing factor such that 

y, &{A f , Z f ) = cr(AA) (7.28) 

Zf 
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Figure 7.1. Nucleon removal number as function of impact parameter in carbon-copper collisions. 

The Rudstam formula for a(AA) was not used because the Ad dependence is too simple and 
breaks down for heavy targets (Townsend et al. 1984; Townsend, Wilson, and Norbury 1985). 

The charge of the removed nucleons A Z is calculated according to charge conservation 

Z P = Z F + AZ (7.29) 


and is divided among the nucleons and alpha particles according to the following rules. The 
abraded nucleons arc those removed from that portion of projectile in the overlap region with 
the target. Therefore, the abraded nucleon charge is assumed to be proportional to the charge 
fraction of the projectile nucleus as 


Z;ihr 


Zp Aabr 

A P 


(7.30) 


This, of course, ignores the charge separation caused by the giant dipole resonance model of 
Morrissey et al. (1978). The charge release in the ablation is then given as 


^abl 


— AZ — Z abr 


(7.31) 


which simply conserves the remaining charge. 

The alpha particle is known to be unusually tightly bound in comparison with other nucleon 
arrangements. Because of this usually tight binding of the alpha particle, the helium production 
is maximized in the ablation process 


N a = int 



(7.32) 


where int(x) denotes the integer part of x. The number of protons produced is given by charge 
conservation as 

N P — AZ — 2N a (7.33) 

Similarly, neutral conservation requires the number of neutrons produced to be 


Nn — AA — Np 4A q 


(7.34) 


The fragments with masses of 2 and 3 are ignored. 
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The calculation is performed for A A = 1 to A.4 = Ap — 1, for which the cross section 
associated with A A > Ap - 0.5 is missed. These are, of course, the central collisions for which 
the projectile is assumed to disintegrate into single nucleons if R P < R T as 


N P = Z P (7.35) 

N n = A P - Z P (7.36) 


and is ignored otherwise. The energetic target fragments are being ignored as well as the mesonic 
components. The peripheral collisions with A .4 < 0.5 are also missing. Most important in these 
near collisions is the coulomb dissociation process studied by Norbury and Townsend (1986). 

Only the nuclear radius for use in the model is yet undefined. Equation (7.1) is an accurate 
representation of the high-energy cross sections. The choice of nuclear radius as 


R = l.26A 1 / 3 (7.37) 

is consistent with equation (7.1) for cq = 1.26 fm when the peripheral collisions (A/1 < 0.5) are 
taken into account. This completes the description of the basic fragmentation model in present 
use. 


In the present evaluation, we look only to elemental fragmentation cross sections for which 
most of the experimental data have been obtained. This is also motivated by the crudeness of 
the present model which is not expected to be completely accurate. Even so, the quality of the 
experimental data base is uncertain with experiments of different groups differing by a factor 
of 2, in general, and even more for specific isotopes. 

The first comparison is with the experiments of Heckman (1975) with 12 C ion beams at 
1.05 GeV/amu on the series of targets extending from hydrogen to lead as shown in table 7.1. 
The present calculations are shown as values in parentheses. The calculated values for hydrogen 
targets are those of Rudstam (1966). Note that all values are within 20 percent of the 
experiments with few exceptions (namely, fragments from hydrogen targets and the neutron 
removal cross section in copper and lead targets). 

The charge removal cross sections for several projectiles on carbon targets are given in 
table 7.2. The agreement between the present model and the Lawrence Berkeley Laboratory 
groups (Heckman 1975; Westfall et al. 1979) is quite good. Our results tend to be low compared 
with the experiments of Webber et al. (1983a and 1983b) and Guerreau et al. (1983). The model 
can be adjusted once experimental differences are resolved. 


The elemental fragmentation cross section of iron projectiles on several targets is shown in 
table 7.3. Again, reasonable agreement is found generally throughout the table with a few 
examples of relatively large errors. The quantities in brackets at the bottom of the table are the 
coulomb dissociation cross sections for forming manganese. These are to be added to the nuclear 
fragmentation cross sections for manganese in parentheses before comparing with experimental 
values. 
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Tabic 7.1. Fragmentation Cross Sections of Carbon Beams 
at 1.05 GeV/amu in Various Targets 

[Quantities in parentheses are present theory] 



Carbon cross section,® mb, in target of— 


H 

Be 

c 

Cu 

Pb 

Fragment 

(b) 





Li 

Be 

B 

C 

23 ± 2 (34) 
17 ± 1 (22) 
50 ± 1 (42) 
28 ± 3 (10) 

51 ± 2 (54) 
35 ± 1 (32) 
81 ± 4 (86) 
49 ± 3 (39) 

52 ± 3 (61) 

35 ± 1 (33) 
78 ± 3 (100) 
50 ± 4 (44) 

71 ± 5 (81) 

47 ± 2 (48) 
119 ± 8 (138) 
86 ± 8 (57) 

103 ± 14 (113) 
71 ± 6 (63) 
203 ± 32 (185) 
139 ± 22 (79) 


n Heckman 1975. 

^Values in parentheses in this column are those of modified Rudstam (1966). 


Table 7.2. Charge Removal Cross Sections of Various Projectiles on Carbon Targets 


Quantities in parentheses arc present theory; 
number in brackets is energy in GeV/amu . 


AZ 

Charge removal cross section, mb, of projectile of 


C [2.1] 
(«) 

0 [2.1] 
(“) 

O [0.9] 

(b) 

Ne [0.47] 

(b) 

Ar [0.21] 

(c) 

Ke [1 .88] 

(d) 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

50 ± 4 (40) 
78 ± 3 (100) 
35 ± 1 (33) 
52 ± 2 (61) 

45 ±2 (45) 
105 ± 4 (101) 
116 ±6 (93) 
50 ± 2 (65) 
36 ± 1 (24) 
65 ±3 (47) 


(40) 

(132) 1 

--- (151) 

154 ± 26 (85) 

122 ± 16 (72) 

144 ± 19 (64) 

81 ± 15 (59) 

112 ± 15 (51) 

90 ± 3 (50) 

92 ± 13 (44) 

65 ± 11 (42) 

83 ± 13 (37) 
(35) 

(64) 

181 ± 27 (157) 
124 ± 13 (110) 
100 ±11 (87) 
87 ± 11 (76) 

54 ± 9 (62) 
78 ± 11 (67) 

52 ± 7 (57) 

55 ± 9 (52) 

53 ± 7 (49) 

54 ± 10 (45) 
59 ± 10 (42) 
57 ± 10 (39) 
83 ±11 (36) 
(35) 

176 ± 5 
104 ± 5 
55 ± 3 
27 ± 2 

129 ± 3 (90) 
214 ± 3 (98) 
155 ± 3 (75) 
140 ± 3 (65) 
74 ± 2 (54) 
33 ± 1 (19) 
(40) 



a Heckman 1975. 

^ Webber et al. 1983a and 1983b. 
c Guerreau et al. 1983. 

^Westfall et al. 1979. 


Comparing the model cross sections with the experimental data set reveals that 92 percent of 
the calculated cross sections are within 50 percent of the measured values. If we reduce the error 
band to 30 percent, we will find 81 percent of the cross sections are in agreement to within this 
level. Among the least accurate are the ion on hydrogen target data which again is Rudstam s 
theory and the cross sections of Webber et al. Note that our model agrees with experiments 
to the extent that the experimentalists agree among themselves for the same projectile-target 
combinations. From this point of view, little progress can be made in improving the model 
until the experimental situation is clarified. The model of Silberberg, Tsao, and Shapiro (19 ), 

which includes many corrections to Rudstam’s formulas, is preferred for hydrogen targets. 
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The semicmpirical model for argon fragmentation on carbon is shown with the quantum 
mechanical optical model calculation (Townsend et al. 1986a) in figure 7.2. Also shown are 
experimental data of Viyogi et al. (1979). Reasonable agreement is seen between the two 
models except for neutron removal where there are no data yet to resolve the difference (Wilson, 
Townsend, and Badavi 1987b). 
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• Wilson, Townsend, and Badavi et al. 1987b 

1 Viyogi et al. 1979 
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Figure 7.2. Representative argon-carbon fragmentation cross sections. 


8. Space Radiation Exposure Issues 

The differences in radiosensitivity of various tissues within an individual as well as individual 
differences are generally assumed to result from repair mechanisms (Curtis 1986; Fritz-Niggli 
1988). The work of Swenberg, Holwitt, and Speicher (1990) suggests these differences may result 
from the structural state of the DNA as well. Repair also affects the dose response for protracted 
exposure. Current radiation protection guidelines use quality factors that are independent of 
dose rate (no time modifying factors), which may be of unusual importance in the small dose 
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rate exposures often experienced in space (NCRP 98 (Anon. 1989)). Clearly, well understood 
dose rate dependent models are needed (Curtis 1986; Anon. 1989). Furthermore, exposure 
received on a mission to the Moon or Mars will involve heavy ion exposure for which many issues 
concerning appropriate relative biological effectiveness (RBE) factors (hence, quality factors) are 
yet unknown. The accumulated levels of heavy ion exposure will be large and unprecedented in 
human experience (Nachtwey and Yang 1991). Although these issues may be studied in ground- 
based exposures with model biological systems, extrapolating to human exposure is difficult at 
best and space stress factors such as microgravity are unknown possible modifying factors in 
radiobiological response. 

In addition to radiobiological response issues is the need to evaluate dose nonuniformity 
caused by body self-shielding (Khandelwal and Wilson 1974) and dose gradients within the 
shielding structure. For example, it is well known that tumor prevalence in the female breast is 
site specific even for relative uniform exposure (NCRP 85 (Anon. 1986)). We are led to believe 
that the exposure of only sensitive sites may be effective in tumor formation. Conversely, 
exposure of insensitive sites is assumed to be noneffective, and nonuniformity of exposure is a 
critical issue. In this section, we make some preliminary assessments concerning these issues 
and examine a limited number of shielding strategies to mitigate these radiation effects. 

8.1. Galactic Cosmic Ray Exposure 

The incident galactic cosmic ray spectrum (Adams, Silbcrberg, and Tsao 1981; Adams 1987) 
for free space is propagated through the target material by using the accurate analytical- 
numerical solutions to the transport equation described in chapter 10 of Wilson et al. 1991. 
These solution methods have been verified (to within 2 percent accuracy) by comparison with 
exact, analytical benchmark solutions to the ion transport equation (Wilson and Townsend 
1988; Wilson et al. 1988). 

These transport calculations include 

1. Linear energy transfer (LET) dependent quality factors from ICRP 26 (Anon. 1977) 

2. Dose contributions from propagating neutrons, protons, alpha particles, and heavy ions 
(high energy, high charge (HZE) particles) 

3. Dose contributions resulting from target nuclear fragments, primary particles, and their 
secondaries 

4. Dose contributions due to nuclear recoil in tissue 

Major shortcomings of the calculations are as follows: 

1. Except for tissue targets, contributions of fragments with masses of 2 and 3 are neglected 

2. All secondary particles from HZE interactions are presently assumed to be produced with 
a velocity equal to that of the incident particle; this is conservative for neutrons produced 
in HZE particle fragmentations 

3. Meson contributions to the propagating radiation fields are neglected 

4. Nucleus- nucleus cross sections are not fully energy dependent (nucleon-nucleus cross 
sections are fully energy dependent) 

These items are not conservative and probably alone result in a 15- to 30-percent underesti- 
mate of the exposure. As discussed by Townsend, Wilson, and Nealy (1989) and Townsend and 
Wilson (1988), the main sources of uncertainty are the input nuclear fragmentation model and 
the incident galactic cosmic ray (GCR) spectrum. Taken together, they could easily impose an 
uncertainty factor of 2 or more in the exposure predictions. 
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8.2. Results 


In the present results, we use the ICRP 26 (Anon. 1977) quality factors which are currently 
in force within the U.S. space program. Figure 8.1 displays dose equivalent (in units of sieverts 
per year) as a function of water shield thickness (in units of areal density, grams per centimeter 2 , 
or thickness, centimeters). Curves are displayed for solar minimum and solar maximum periods. 
For all thicknesses considered, the dose and dose equivalent during solar maximum are less than 
half the dose equivalent during solar minimum, at least according to the current estimates using 
the CREME environmental model of Adams (1987). Figure 8.2 displays results for dose and 
dose equivalent behind an aluminum shield. Also shown are measurements using the argon 
ion chambers at two shielded locations (Kovalev, Muratova, and Petrov 1989). The results 
for the location at 1 g/cm 2 are the most clear by experimental design. The mass distribution 
for the deeply shielded counter was poorly defined (Kovalev, Muratova, and Petrov 1989); this 
uncertainty is denoted in the figure by using parentheses around the data points. The solar 
maximum model predicted by CREME is clearly an underestimate. The solar minimum model 
appears in reasonable agreement with the Prognoz spacecraft data. Therefore, we will restrict 
the present analysis to solar minimum periods, which are the most limiting for GCR exposures. 
This does not imply, however, that exposures during solar maximum periods are not important. 
On the contrary, the cumulative exposures resulting from combined GCR and increased solar 
flare activity during solar maximum could potentially be significant (Nealy et al. 1990). During 
solar minimum periods, the estimated unshielded dose equivalent of 1.2 Sv does not exceed the 
exposure limits for either the skin or the ocular lens (which are 3 Sv for skin and 2 Sv for ocular 
lens). The dose equivalent at a depth of 5 cm, which yields an estimate of the exposure to 
the unshielded blood-forming organs (BFO), is 0.61 Sv, which exceeds the limit of 0.5 Sv by 
22 percent. To reduce this estimated exposure below 0.5 Sv requires approximately 3.5 g/cm 2 
(3.5 cm) of water shielding in addition to the body self-shielding of 5 g/cm 2 (5 cm). 



Dose equivalent 
Theory 

£ Prognoz (Kovalev, Muratova* and Petrov 1989) 


Dose 


Theory 

i Prognoz (Kovalev, Muratova* and Petrov 1989) 



Figure 8.1. Dose equivalent resulting from galac- 
tic cosmic rays as function of water shield 
thickness. 


Figure 8.2. Deep space exposure behind alu- 
minum shield. Parentheses denote depth in 
interior of Prognoz spacecraft. 


For relative comparison purposes, the BFO dose equivalent as a function of shield thickness 
(areal density) is plotted in figure 8.3 for three materials (aluminum, water, and liquid hydrogen). 
Shielding effectiveness per unit mass increases as the composition of the shield changes from 
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heavier to lighter mass elements. For liquid hydrogen, an added advantage is the reduced neutron 
fluence caused by the absence of neutrons in the target composition and by the lack of target 
fragment contributions because of the elementary nature of hydrogen. From these results, for 
an allowed BFO exposure of 0.25 Sv/year, which corresponds to an uncertainty factor of 2 in a 
0.5 Sv/year estimate, the mass ratios for the shielding are approximately 1:5:11 for LH2:H20:A1. 
Obviously, for GCR shielding, the materials of choice are those composed of low atomic mass 
number constituents with significant hydrogen content. 



Figure 8.3. Dose equivalent in BFO as function of shield type and thickness. 


Although the calculations are useful for estimating relative shield effectiveness to compare 
different materials, quantitatively the calculations should be considered as preliminary estimates 
of actual shield mass requirements. Aside from the previously mentioned shortcomings related 
to neglecting meson production and target fragment contributions from interactions of HZE 
particles and the target medium, figure 8.3 shows that the dose equivalent is a slowly decreasing 
function of shield thickness. This is a result of secondary particle production processes whereby 
the heavier GCR nuclei are broken into nucleons and lighter nuclear fragments by nuclear and 
coulombic interactions with the shield material. This slow decrease in dose equivalent with 
increasing shield thickness means that relatively small uncertainties in predicted doses arising 
from nuclear fragmentation model inaccuracies may yield large uncertainties in estimated shield 
thicknesses. A preliminary analysis of the nonlinear relationship between exposure uncertainty 
was presented by Townsend, Wilson, and Nealy (1989). The most startling finding was that an 
uncertainty factor of 2 in exposure amplified into an order of magnitude uncertainty in shield 
mass requirements. 

8.3. Biological Effect Modifications 

The rising RBE at low GCR dose results from the multitarget assumption in Katz theory 
leading to the sigmoid behavior in the survival curve of low-LET radiation as opposed to the 
exponential relationship for high-LET radiation (Cucinotta et al. 1991). The transition from 
sigmoid to exponential behavior is observed by Todd and Tobias (1974) to occur at 150 to 
200 keV/pm for mammalian cells. Many also believe that the sigmoid behavior is related to 
repair mechanisms. This view is promoted by single exposure and split exposure experiments 
with a delay of 2.5 or 23 hours between fractions using V79 hamster cells as shown in figure 8.4 
(Elkind and Sutton 1960). Repair is indicated by the sigmoid response of the second exposure 
after either the 2.5-hour repair period or especially the 23-hour repair period. Obviously the 
RBE based on such a photon exposure protocol depends on the history of the radiation induced 
damage. Similar survival studies with confluent C3H10T1/2 mouse cells (Gl) indicate no repair 
for this end point for high-LET radiations (Yang et al. 1986). As a result of operative repair 
mechanisms (sparing) for low-LET exposure (fig. 8.5) and the lack of repair for high-LET 
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exposure (fig. 8.6), the corresponding RBE is dose rate dependent (Yang et al. 1986) as shown 
in figure 8.7. Also shown in figure 8.7 are the RBE values for neoplastic transformations. 
(Note, great liberty has been taken in connecting the data points.) The increase in RBE at 
low dose rate is in part indicative of repair of the damage for low-LET radiation (fig. 8.8) but 
additional enhancement of high-LET exposure at low dose rate (presumably some misrepair 
mechanism) also contributes for cell transformations (fig. 8.9). If misrepair/ repair plays a role 
then this should be observed in the delayed plating experiments of Yang et al. (1989) as shown 
in figure 8.10. Instead the delayed plating experiments show no transformation misrepair but 
repair appears in cell survival data in distinction to the earlier low dose rat(> experiences with 
the same cell system. 



Figure 8.4. Fractional survival of cultured Chi- 
nese hamster cells for single exposure and ex- 
posure in two fractions (Elkind and Sutton 
1960). 



Figure 8.5. Survival fraction of confluent mouse 
cells at two dose rates displaying sparing at 
low dose rate (Yang et ah 1986). 


Similar dose rate enhancement effects are observed in asynchronous cell cultures by Hill et al. 
(1982 and 1985, fig. 8.11) and whole animal exposures as observed by ThontSon et al. (1981a 
and 1981b), Thomson, Williamson, and Grahn (1983, 1985a, 1985b, and 1986), and Thomson 
and Grahn (1988 and 1989) (fig. 8.12). These effects are considered the result of cell cycle 
phenomena (Rossi and Kellerer 1986; Brenner and Hall 1990). The basic model assumes that 
some phases of the cell cycle are more affected by radiation exposure. This is clearly seen in 
the cell synchronous experiments of Terasima and Tolmach (1963) shown in figure 8.13. The 
model of dose rate enhancement assumes only one cell phase is effective in injury of only that 
fraction in the sensitive phase. At a later time, a different fraction of cells is in the appropriate 
phase providing two exposed groups of cells and an apparent enhancement. Such a model 
was exploited in the work of Brenner and Hall (1990). This explanation fails to explain the 
enhancement effects observed by Yang et al. (1989) in cell transformation in stationary phase 
(Gl) confluent C3H10T1/2. 
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Figure 8 . 6 . Dose rate effects on confluent mouse 
cell survival for high-LET exposure (Yang 
et al. 1986). 
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Figure 8 . 8 . Repair processes for confluent mouse 
cell cultures exposed to 7 -rays at low LET. 
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Figure 8.7. RBE as function of dose rate. The 
curves are to guide eye. 



Dose, Gy 


Figure 8.9. Cell transformation rate enhanced 
at low dose rate for high-LET exposure with 
possible misrepair mechanism indicated. 





Figure 8.H. Transformation frequencies in C3H10T1/2 cells. 
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Figure 8.12. Life shortening in mice after single, 
24, or 60 fractions of neutrons and after sin- 
gle, 24, or 60 fractions of y-rays. Curves are 
fit to data of Thomson et al. 



Figure 8.13. Survival of synchronized human cells 
exposed to 500 rads of 200 kV X-rays. 


Clearly the risk to long-term GCR/SCR exposure will be difficult to evaluate because of 
the low dose rate, fractionated components, and the complex mixture of low- and high-LET 
radiations in space. Operative repair and cycle enhanced effects will require at least an intimate 
knowledge of the LET distributions at affected tissues or possibly more comprehensive track 
structure data. 


8.4. Nuclear Models, Materials, and LET Spectra 


As is clear from section 8.3, the distribution of exposure component over LET is a primary 
indicator of biological response. For example, low-LET components allow certain biological 
repairs at low dose rates and a low RBE value, whereas high-LET components can show increased 
biological effects at low dose rates and generally high RBE values. There is clear evidence that 
the relative contributions to exposure from various LET components can be altered through 
the choice of shield material. The transmitted LET spectrum for an aluminum shield is shown 
in figure 8.14, and the transmitted LET spectrum for a liquid hydrogen shield is shown in 
figure 8.15. Note that the LET spectra in figures 8.14 and 8.15 are for LET in the respective 
shield material. Although a rather large shift in LET can be accomplished by choice of shield 
composition, an exact evaluation must await improved nuclear fragmentation cross sections 
since uncertainty in cross sections cause LET shifts of the same order of magnitude. These 
shifts can be seen when the LET spectra in the Earth’s atmosphere obtained by using the 
Bowman, Swiatecki, and Tsang (1973) fragmentation model shown in figure 8.16 are compared 
with the spectra obtained by using the Langley fragmentation model shown in figure 8.17. 
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Figure 8.14. GCR integral LET spectra in alu- 
minum for 30° orbit at altitude of 400 km. 


Figure 8.15. GCR integral LET spectra in hydro- 
gen for 30° orbit at altitude of 400 km. 




Figure 8.16. GCR integral LET spectra in 
Earth’s atmosphere for fragmentation pa- 
rameters of Bowman, Swiatecki, and Tsang 
(1973) model. 


Figure 8.17. GCR integral LET spectra in 
Earth’s atmosphere for fragmentation param- 
eters of Langley Research Center model. 


8.5. Human Geometry Factors 

The significance of improving the accuracy (Shinn, Wilson, and Nealy 1990) for predicting 
the dose and dose equivalent that astronauts will incur during future NASA missions has 
been demonstrated in several studies (Townsend, Wilson, and Nealy 1988; Nealy, Wilson, and 
Townsend 1989). For example, Townsend, Wilson, and Nealy (1988) indicate that an increase of 
20 percent in predicted BFO dose equivalent due to GCR s equates to a tripling of the required 
shield mass from 5 to 16 g/cm 2 of water to meet the recommended annual BFO limit of 0.5 Sv. 
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Large uncertainties are associated with the current dose estimate analysis, and every possible 
effort is needed to improve the accuracy to accomplish these missions in the most economical 
way without compromising the well-being of the astronauts. 

One of the customary estimation practices (Beck, Stokes, and Lushbaugh 1972) that has 
been considered fairly reliable in the past is the use of the equivalent sphere model to obtain 
dose or dose equivalent to BFO. Langley and Billings (1972b) examined the feasibility of using 
a set of dosimetry spheres to monitor real-time organ doses received by astronauts under 
various space radiation and vehicle conditions. They made comparisons between the doses 
calculated for the spheres and the detailed body geometry under a range of solar proton energy 
spectrum characteristics and also under various vehicle radiation shielding thicknesses. The 
spectra were characterized by an assumed form described by Webber (1963). The optimal 
radii were determined for those spheres with the corresponding correlation constants that best 
represented the averaged organ doses under those assumed conditions. Although a moderate 
error of 18 percent for the correlation was found, one might question if the accuracy will hold 
under less idealized particle spectral conditions. This question was considered by Shinn, Wilson, 
and Nealy (1990) and we discuss those results herein. 

The calculation made in a separate study (Simonsen et al. 1990) for the radiation transport 
through the atmosphere of Mars for the three largest solar flares observed in the last half 
century is extended here to include detailed BFO geometry. Comparisons are made for the dose 
equivalent to the various distributed BFO with the reported values based on the equivalent 
sphere model. 

The Langley Research Center nucleon transport code BRYNTRN (Wilson et al. 1989b) was 
used by Simonsen et al. (1990) to obtain dose and dose equivalent on the surface of Mars caused 
by large solar flares. The transport code was based on the straight ahead approximation, which 
reduces consideration to one-dimensional transport; the merits of this approximation have been 
discussed elsewhere (Alsmiller et al. 1965; Alsmiller, Irving, and Moran 1968). An asymptotic 
expansion for the solution to the transport equation in two dimensions, subject to boundary 
conditions given for an arbitrary convex region, was derived by Wilson and Khandelwal (1974). 
The first term of the expansion was found to be an accurate approximation of the dose and for 
the case of an isotropic proton fluence spectrum is given by 


with 


D(x) = 47r 


roc r o 

/ 4>{E) / 

J 0 Jo 


R(t , E) f x (t) dt dE 



dt = 1 


( 8 . 1 ) 


( 8 . 2 ) 


where R(t , E ) is the fluence-to-dose conversion factor at the depth t for normal incidence protons 
on a slab and f x (t) is the areal density distribution function for the point x. The quantity 
fx(t) dt is the fraction of the solid angle for which the distance to the surface from the point x 
lies between t and t + dt. 

To simplify the computational task (that is, without making any change to the BRYNTRN 
code), equation (8.1) is rewritten as 
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where D x (t) is the dose (or dose equivalent) at depth t for normal incidence protons on a slab 
of tissue. With the areal density distribution function for BFO given by the detailed geometry 
work described by Langley and Billings (1972a) and Billings and Yucker (1973), equation (8.3) 
can be calculated. 

8.6. Results 

The three solar flare spectra used for this study are those of February 1956, November 1960, 
and August 1972 events, whereas Langley and Billings used a Webber (1963) form of integral 
spectra given by the inverse exponential of proton magnetic rigidity with a range of rigidity 
parameter P D from 50 to 200 MV. Figure 8.18 shows these three flare spectra and the best fit 
to the earlier two events with the Webber form. The actual spectra, especially the high-energy 
range of the February 1956 event (Foelsche et al. 1974), are different from the analytical form of 
Webber. The actual spectrum (Wilson and Derm 1976) for the August 1972 event is accurately 
approximated by the Webber form. 



(a) February 1956 event. (b) November 1960 event. 



(c) August 1972 event. Webber spectrum fits actual event spectrum. 

Figure 8.18. Fluence spectra for three major solar particle events. 

The average dose equivalents at the surface of Mars caused by these three solar flare events are 
shown in figure 8.19 as a function of slab (water) thickness for the low-density Mars atmosphere 
model (16 g/cm 2 C0 2 vertically) used in Simonsen et al. (1990). These average dose equivalent 
values are obtained by summing the directional (anisotropic) dose equivalent over the solid angle 
and are used as D x (t) in this section. The calculated results from equation (8.3) are presented 
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in table 8.1 for the five distributed compartments of the blood-forming organs. Also shown for 
comparison are the average BFO and 5-cm (water) depth dose equivalents. 



Figure 8.19. Dose equivalent at Mars surface as function of slab (water) thickness for low-density Mars atmosphere 

model. 

It is customary (Space Science Board 1970; Beck, Stokes, and Lushbaugh 1972) to represent 
the average BFO exposure (dose or dose equivalent) with the 5-cm sphere based on the 
recommendation of the Space Science Board (1970). Conversely, the average BFO dose was 
found to be approximately half the 5-cm sphere dose in several analytical findings, such as the 
one from Langley and Billings (1972b). For the August 1972 event, the average BFO value 
for the detailed geometry (table 8.1) is fairly close (within 10 percent) to one half the value 
for a 5-cm sphere. However, the differences are larger for the other two flares, with 30 and 
41 percent for November 1960 and February 1956 spectra, respectively. This wide discrepancy 
among these three events probably occurs because the two earlier flares contain more penetrating 
high-energy protons (fig. 8.18) and the actual spectra do not conform to the simple analytical 
form that Langley and Billings (1972b) used. We further note that the 5-cm sphere dose is 
conservative for these three events. Also, the 5-cm sphere dose is conservative as was found for 
GCR exposure by Townsend, Shinn, and Wilson (1991). 


Table 8.1. BFO Dose Equivalent at Mars Surface for Low-Density Atmosphere Model 


Solar 

flare event 



Dose equivalent, cSv, for 

— 


Arms 

Legs 

High trunk 

Low trunk 

Skull 

Average 
BFO value 

5-cm sphere 

Feb. 1956 

8.74 

8.60 

8.32 

7.98 

8.91 

8.45 

9.94 

Nov. 1960 

5.66 

5.34 

4.95 

4.32 

5.75 

5.21 

7.31 

Aug. 1972 

3.20 

2.73 

2,42 

1.76 

3.09 

2.56 

4.61 
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9. Concluding Remarks 

Great progress has been made in developing codes for space shielding for future NASA 
programs. It is likewise clear that major uncertainties remain in the environmental model, the 
nuclear cross sections, and the methods for estimating biological risk. These uncertainties have 
an important impact on shield design and mission cost. Substantial work remains before all 
these issues can be resolved. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
November 8, 1991 
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